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We investigate the phase structure of two- component Bose-Einstein condensates (BECs) with 
repulsive intra- and interspecies interactions in the presence of respective vortex-antivortex super- 
position states (VAVSS). We show that different winding numbers of vortex and antivortex and 
different intra- and interspecies interaction strengths may lead to different phase configurations, 
such as fully separated phases, inlaid separated phases, asymmetric separated phase, and partially 
mixed phases, where the density profile of each component displays a petal-like (or modulated petal- 
like) structure. A phase diagram is given for the case of equal unit winding numbers of the vortex and 
antivortex in respective components, and it is shown that conventional criterion for phase separation 
of two-component BECs is not applicable for the present system due to the VAVSS. In addition, our 
nonlinear stability analysis indicates that the typical phase structures of two-component BECs with 
VAVSS allow to be detected in experiments. Moreover, for the case of unequal winding numbers of 
the vortex and antivortex in respective components, we find that each component in any of the pos- 
sible phase structures is in a cluster state of vortices and ant ivort ices, where the topological defects 
appear in the form of singly quantized visible vortex, or hidden vortex, or ghost vortex, depending 
on the specific parameters of the system. Finally, a general rule between the vortex-antivortex 
cluster state and the winding numbers of vortex and antivortex is revealed. 

PACS numbers: 03.75.Mn, 03.75.Lm, 05.30.Jp, 67.85.Fg 



I. INTRODUCTION 



Multicomponent systems are of fundamental impor- 
tance in many fields of physics on account of their univer- 
sality and complexity. In particular, great attention has 
been recently paid to a mixture of two-component Bose- 
Einstein condensates (BECs) in cold atom physics due to 
ultra high purity, excellent theoretical description, exper- 
imental accessibility and precise controllability of a con- 
densate system [1, 2]. Two-component BECs can consist 
of two different alkalis, or different isotopes, or different 
hyperfine states of the same alkali atom. Such a mix- 
ture provides an attractive and versatile testing ground 
to study the intriguing properties of macroscopic quan- 
tum many-body systems inaccessible in single-component 
BECs. In fact, many interesting phenomena have been 
predicted theoretically and some observed experimen- 
tally in two-species BECs, ranging from the topologi- 
cal structures of the ground and excited states [3l-[l^, 
symmetry-breaking transition p^], quantum turbulence 
[16], pattern formation vortex-bright-solitons [18| 

and skyrmions [l9|, entangled states [20], to the stripe 
phase induced by spin-orbit coupling etc. Among 
these phenomena, of crucial significance is the phase 
structure of the system [sl-fl^. For instance, a separated 
phase of a two-component BEC is the premise of study- 
ing the Kelvin-Helmholtz instability [22] and Rayleigh- 



Taylor instability [23|. 

To the best of our knowledge, most of theoretical litera- 
tures use the Thomas-Fermi (TF) approximation (a semi- 
analytical method) or a variational approach to study the 
ground state properties of two-component BECs, where 
spherically symmetric mixed phases, spherically symmet- 
ric separated phases, and asymmetric side-by-side sep- 
arated phases are predicted. For the former case, the 
kinetic energy is completely neglected. However, when 
the contribution of kinetic energy to the total energy is 
considerable, the TF approximation will fail ^ and thus 
can not be relied upon to determine the phase struc- 
ture of a two-component BEC. The effect of kinetic en- 
ergy on the phase transition of a two-component BEC 
in an infinitely deep square well potential has been re- 
cently demonstrated in Ref. [ll|. As for the variational 
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approach, it is usually limited to few simple quantum 
states and certain external potentials. In this context, it 
is necessary to resort to numerical methods for obtaining 
reliable phase structures. 

In this paper, we investigate the exact two-dimensional 
(2D) phase structures of two-component BECs with re- 
spective vortex-antivortex superposition states (VAVSS). 
The VAVSS in BECs is a quite interesting research object 
as it can exhibit rich physical properties [24-28], such as 
pec uliar petal-like structure [M^^S] and unique dynamics 
|28j . Particularly, the equilibrium properties of a quasi- 
2D degenerate boson-fermion mixture (DBFM) with a 
bosonic VAVSS is recently studied by using a quantum- 
hydrodynamic model [29j. It is shown that the VAVSS 
greatly influences the equilibrium state and stability of a 



2 



DBFM, where a special intermittency phenonienon exists 
in the two stability curves of the DBFM with a bosonic 
VAVSS [m and the one with a bosonic vortex [30|. On 
the other hand, the VAVSS may have many potential ap- 
plications in quantum information, quantum communica- 
tion, and inertial sensing 0, HH, 113 • Experimentally, the 
creation of VAVSS in BECs has been reported by differ- 
ent groups [31-34]. An interesting question is to ask how 
the phase structure of a two-component BEC is modified 
by the respective VAVSS combining with intra- and in- 
terspecies interactions. We show that, depending on the 
relative interaction strengths and the winding numbers 
of vortex and antivortex, the two-component BECs with 
respective VAVSS can display rich phase structures not 
met in other systems, such as fully separated phases, in- 
laid separated phases, and asymmetric separated phase 
with (deformed) petal-like component density profiles. 
We give a phase diagram for the case of equal unit wind- 
ing numbers of vortex and antivortex in respective com- 
ponents. The typical phase structures are long-lived ac- 
cording to our nonlinear stability analysis. For the case of 
unequal winding numbers of vortex and antivortex, each 
component in any of the possible phases is in a vortex- 
antivortex cluster state, where the phase defects emerge 
in the form of visible vortex, or hidden vortex, or ghost 
vortex. Furthermore, we reveal and discuss the general 
relation between the vortex- antivortex cluster state and 
the winding numbers of vortex and antivortex. 

The paper is organized as follows. In Sec. II, we de- 
scribe the model for two-component BECs with respec- 
tive VAVSS. In Sec. Ill, we study the phase structure of 
two-component BECs with respective VAVSS with equal 
unit winding numbers of the vortex and antivortex. A 
phase diagram is presented, and the stability of phase 
structures is analyzed. In Sec. IV, we discuss the phase 
structure of the system with unequal winding numbers 
of the vortex and antivortex in respective components. 
The conclusion is outlined in the last section. 



II. MODEL 

We consider a 2D system of trapped two-component 
BECs with tight confinement in the z direction. In the 
zero-temperature hmit, the dynamics of the system is de- 
scribed by the coupled Gross-Pitaevskii (CP) equations 



-vi + y(r)+5ii^/;ir+5iV2rj ^/'i, (1) 

-Vi + y(r) + 52|^/'2|' + <7|Vi|'l ^/'2, (2) 



where we have adopted plane-polar coordinates for space 
(r,0), Vi = i|:(r|:) + and V{t) = r^A is the 

trapping potential. Here we have assumed that the two 
species have the same atomic mass m and undergo the 
same trapping frequency uj^ in the x-y plane. In this pa- 
per, length, ti me, energy , and angular momentum are in 
units of do = \/h/2muj^^ fvuj^^ and respectively. 



^1,^2, and g denote dimensionless intra- and intercompo- 
nent coupling strengths which are proportional to the 
corresponding 5- wave scattering lengths aj^aij {j = 1,2) 
between intra- and intercomponent atoms and the atom 
numbers Ni^2 in species 1 and 2 [1, 2, 28^ _29]- The two 
wave functions are normalized as // |'0i,2| rdrdO = 1. 
We seek dimensionless solutions of stationary VAVSS 

^j{r,0,t) = ^j{r,0)e-'^^', j = l,2, (3) 

where /ij is the chemical potential, and cpj is given by 



J = 1,2. (4) 



Here fjk{r) = Aje-"^"/^-" {r/ajY^" and fjp{r) = 

Aje~^ /^^j {r/ajY^p are the amplitudes of the vortex and 
antivortex components with Aj being the normalization 
constant and aj the width of condensate j. Ijk and Ijp 
are the winding numbers of the vortex and antivortex, 
and the real constants aj , Pj show the proportion of the 
vortex and antivortex with + = 1. The relative 
phase Sj just causes offset of the density profile by an 
angle Sj/{ljk + Ijp)- Note that in Eq. (j4|) we have taken 
into account the infiuence of winding number on the con- 
densate size in terms of and a - . 

Substituting Eq. ([3]) into Eqs. ([1]) and ([2]), we obtain 
the coupled equations for (^i and (/?2, 



rdr dr ^2 ikj 
+%^^ + {V + gM^+g\^2?)^i, (5) 



I d dip2 ^ f2ph 

r dr dr 

-%^2 + (V^ + 52|^2|' + 5|¥'i|')^2, (6) 



where the second and third terms associated with wind- 
ing numbers in the right-hand sides of Eqs. (j5j) and (j6j) 
are resulted from the presence of VAVSS in each compo- 
nent. Consequently, the chemical potentials /ii,2 read 



Ml 



2(/ifc + l) 



1 



\i^i\^r^drde 



47r(/ip - hk) 



j flA 



rdr 



+ IJiV ^ gil^if ^ g\^2f]\^if rdrdO, (7) 

2{l2k + 1) 1 //, ,2 3. 

2 4 / / 1^2 1 r^drdO 

^2 ^2 J J 

flpl3lrdr 
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[V + g2\^2f + g\^if]\^2f rdrde, (8) 



and they must be determined self-consistently from the 
normalization conditions of (fi and (/?2- 
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From Eqs. (j5j)-(j8j) we can see that when the kinetic 
energy is comparable with the sum of interaction energy 
and potential energy the TF approximation is not suit- 
able to describe the system, especially for cases of VAVSS 
with large winding numbers. In the following, we numer- 
ically solve the coupled equations (j5j) and ([6]). Starting 
with two trial VAVSS with specific ratios of Oj and 
we obtain the exact 2D equilibrium state of the system by 
using the imaginary time propagation method based on 
the Peaceman-Rachford method [35] , which is equivalent 
to the procedure of minimizing the sum of the average 
energy per atom in species 1 and that in species 2, 

E = JJ rdrdO[i^i \(pif ^ \^2f - ^gi 

-^92 1^2^ -gl^lf \^2f]- (9) 

III. STRUCTURE OF TWO-COMPONENT 

BECS WITH RESPECTIVE VAVSS WITH 
EQUAL WINDING NUMBERS OF VORTEX 
AND ANTIVORTEX 

Here we just consider the case of = = 1/2 be- 
cause the petal structure of a VAVSS with unequal su- 
perposition ratio only emerges under the condition of a 
very small particle number or an extremely weak inter- 
atomic interaction which is usually not met in exper- 
iments [29|. For convenience, we introduce two relative 
interaction strengths, R21 = ^2/^1 and R = g/gi^ and as- 
sume the intra- and interspecies interactions to be repul- 
sive. Throughout this paper the relative phase is taken to 
be 5j =0. Figure 1 shows the phase diagram of two BECs 
with VAVSS with lik = hp = ^ and hk = hp = where 
gi = 200. There exist eight possible phases depending on 
the values of R21 and i?, and the typical density profiles 
corresponding to phases I- VIII are displayed in Fig. 2 
and Figs. 3(a)-3(b), respectively. In any of these phases, 
the density profile of each component displays a petal 
structure or a crescent-pair (deformed petal) structure 
or a combined structure of petal and crescent-pair due to 
the presence of respective VAVSS and the competition 
between the intraspecies and the interspecies repulsions. 
In Fig. 1, regions I and IV denote two opposite fully 
separated phases. In the former phase, component 1 is 
completely expelled outside component 2 [see Figs. 2(al) 
and 2(a2)], while in the latter phase the sequence is con- 
verse [Figs. 2(dl) and 2(d2)]. Regions II and HI repre- 
sent two opposite inlaid separated phases, where in the 
former case species 1 lies separately in the interval be- 
tween the inner petal layer and the outer crescent layer 
of species 2, however in the latter case the arrangement 
of the two species is inverted [Figs. 2(bl)-2(c2)]. Regions 
V and VI are two different partially mixed phases, where 
in the first phase species 2 is merged inside species 1 
but it is reversed in the second phase [Figs. 2(el)-2(f2)]. 
Finally, the vertical solid-line section {0 < R < 3.5) of 
R21 = 1 marks a fully mixed phase in which the mixing 
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FIG. 1: (color online) Phase diagram of two- component BECs 
with VAVSS with Ijk Ijp = 1 (j = 1,2), where R21 
92/91 and R — 9/91. The parameters are gi — 200 and 
= /^^ = 1/2. Regions I and IV, II and III, and V and VI 
represent two different fully separated phases, inlaid separated 
phases, and partially miscible phases, respectively. Vertical 
solid-line section VII marks a fully mixed phase, while dashed- 
line section VIII denotes an asymmetric separated phase. 



reaches the maximum and the densities profiles of two 
components are the same [Figs. 3(al) and 3(a2)], while 
the vertical dashed-line section {R > 3.5) signs an asym- 
metric separated phase in which both the density profiles 
are symmetric breaking upon the trap center [Figs. 3(bl) 
and 3(b2)]. 

The ground state structure of two-coniponent BECs 
has been studied by several theoretical [3|-[8|, [lo|, [Hi and 
experimental groups [13, 14]. A conventional criterion for 
phase separation is given by gig2 < which is based 
on the minimization of the total interaction energy [l|, Q . 
Note that in the derivation of this criterion the BECs 
are supposed to be homogeneous and the kinetic energy 
is neglected. However, the criterion will fail when the ki- 
netic energy becomes important in a nonuniform system. 
In fact, the kinetic energy play a vital role in determining 
the configuration of two-component BECs with VAVSS. 
According to the condition 9192 < g^^ for instance, our 
system would be in a separated phase if R21 = 2 and 
R> V2, but only R > 2.58 can the actual phase separa- 
tion occur when R21 = 2. The similar effect arises for any 
fixed value of R21 or R as shown in Fig. 1. Physically, the 
kinetic energy acts against the interspecies interaction. 
The latter is responsible for phase demixing while the 
former tends to expand the BECs and thus favors phase 
mixing. At the same time, the trapping potential tends to 
trap the condensates and hence also sustains phase mix- 
ing. Therefore, phase separation can be suppressed by 
the kinetic energy and external potential in some situa- 
tions even if the condition 9192 < g^ is satisfied, as we saw 
in Figs. 1-3. Recently, the influence of kinetic energy on 
the mixing- demixing transition by changing the confine- 
ment was discussed in Ref. pjj]. The authors introduce 
a parameter r] = J (i7^?/^i?/^2 to characterize the overlap 
between two condensate wave functions. And then the 
system shows phase separation if 77 <C 1. However, nei- 
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FIG. 2: (color online) Density profiles (a)-(f) correspond to 
phases I- VI in Fig.l, respectively. Here 1 and 2 label two 
different components of the system, (a) R21 = 0.2, R = 8,(b) 
i?2i = 0.8, R = 8,(c) R21 = 1.5, i? = 7,(d) R21 = 2.5, = 
5,(e) R21 = 0.6, R = l,and (f) R21 = 1.4, R = 1.6. The other 
parameters are the same as those in Fig.l. The darker color 
area indicates the lower density, x and y are in units of do. 
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FIG. 3: (color online) (a)-(b): Density profiles (al)-(a2) and 
(bl)-(b2) correspond to phases VII and VIII in Fig.l, respec- 
tively, where (a) R21 = 1,R = 2 and (b) R21 = l,i? = 8. 
The other parameters are the same as those in Fig.l. (c)-(d): 
Temporal evolution of density profiles, where the initial state 
is given by Fig. 2(a) with a random perturbation. The time 
is (c) t = 0.4 and (d) t = 4. Here x and y are in units of do, 
and t is units of 1/(jU±. 



ther the criterion gig2 < nor the condition 77 <C 1 can 
distinguish different separated phases as shown in Fig. 2 
and Fig. 3(b). In addition, from Eqs. ([I]), (|2|), (|5]) and 
(|6|), we can see that the two components of the system 
satisfy the exchange symmetry, which indicates that the 
phase diagram will be the same if one exchanges the two 
component wave functions. For fixed values of gi and g, 
nevertheless, two different values of ^2 which are sym- 
metric concerning the line of R21 = 1 will lead to two 
different quantum states. The corresponding symmetric 
one of the interaction set (^1,^2,^) is (^2,^1,^) [si swaps 
^2, and the set (^2,^1,^) is not in the same phase dia- 
gram] rather than the set (^1, 2^1—^2, ^) being symmetric 
about the line of R21 = 1. This point can explain why 
the phase diagram in Fig. 1 is asymmetric with respect 
to the line of R21 = 1. Here gi = 200 is a typical pa- 
rameter value in current experiments. If the value of gi 
does not vary largely, the phase diagram will be basically 
unchanged except for possible little offset of the bound- 
ary Hues between different phases. When the value of gi 
changes largely, there exist similar phase structures and 
the regions corresponding to different phase structures 
will possibly redistribute in the phase diagram. 

The phase structures of two-component BECs with re- 
spective VAVSS are obviously different from those of the 
usual two-component BECs [3-14] by virtue of the petal 
structure of the VAVSS. On the other hand, the equilib- 
rium properties of the present system are also evidently 
different from those of DBFM with a bosonic VAVSS [29] . 
In the latter case, due to the Pauli exclusion principle, 
there is no s-wave interaction between identical fermions 
in the spin polarized state. In addition, for the sepa- 
rated phases in the latter case, the gap region between 
two bosonic petals is always occupied by the Fermi gas. 



while for the fully separated phases and inlaid separated 
phases in the former case there is no particle occupation 
in the gap region of two petals of any component due to 
the interspecies repulsion and the requirement of sym- 
metric petal structure and angular momentum conserva- 
tion. Furthermore, here the inlaid separated phases and 
the asymmetric separated phase have no counterparts in 
DBFM with a bosonic VAVSS. 

The VAVSS is a collective excitation (a metastable 
state) of BECs, which is similar to the case of a sin- 
gle pure vortex state. Although the energy of a VAVSS 
is higher than that of a pure vortex [25], the configura- 
tions of two BECs with VAVSS may be long-lived. To 
verify this point, we perform a nonlinear stability anal- 
ysis by monitoring numerically the evolution of a per- 
turbed stationary solution in Eqs.(j5j) and ([6j). Shown in 
Figs. 3(cl)-3(d2) are the time evolution of density pro- 
files, where the initial state is fully separated phase I [see 
Fig. 2(a)] with a random perturbation. The system will 
collapse when the initially density profiles begin to de- 
form severely. This method can effectively estimate the 
Hfetime of a VAVSS or a phase structure of the system. 
As a matter of fact, the nonlinear stability analysis is 
widely used in the study of solitons [36]. Our simula- 
tion shows that the lifetimes of phases I- VII can reach 
the order of 100 ms and the lifetime of asymmetric sep- 
arated phase VIII is a bit shorter due to the symmetry 
breaking, which allows to be detected in experiments. It 
is deserved to mention that similar phase structures also 
exist in two-component BECs with VAVSS with higher 
winding numbers lik = hp = hk = hp- In that case, the 
stability of phase structures become lower because of the 
higher energy resulted from the higher winding numbers. 
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FIG. 4: (color online) Density profiles of two- component 
BECs with respective VAVSS with lik = /2fc = 1 and 
/ip = l2p = 2. (a) i?2i = 0.2, R = 8,(b) i?2i = 0.9, i? = 8,(c) 
i^2i = 1.1, = 7,(d) R21 = 2.5, = 5,(e) R21 = 0.6, = 
l,and (f) R21 = 1.4, = 1.6. The other parameters are the 
same as those in Fig. 2. Here x and y are in units of do. 
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FIG. 5: (color online) Phase distributions corresponding to 
those of the density profiles of Fig. 5. The value of the phase 
varies continuously from — tt to tt. The dark color area indi- 
cates the lower phase. Here x and y are in units of do. 



IV. STRUCTURE OF TWO-COMPONENT 
BECS WITH RESPECTIVE VAVSS WITH 
UNEQUAL WINDING NUMBERS OF VORTEX 
AND ANTIVORTEX 

In Fig. 4, we present the typical structures of two- 
component BECs with respective VAVSS with Ijk = 1 
and Ijp = 2 (j = 1,2). Depending on the intra- and 
interspecies coupling strengths, the system shows differ- 
ent structures such as fully separated phases, inlaid sepa- 
rated phases, and partially mixed phases, which is similar 
to the case of Ijk = Ijp = 1 to some degrees. The density 
profile of each species in a fully mixed phase is similar to 
Fig. 4(e2) or Fig. 4(fl). Here an interesting characteris- 
tic is that the density profiles, except those of the inner 
species in partially mixed phases and the two species in 
a fully mixed phase, generally form closed (highly mod- 
ulated petal-like) structures. The corresponding phase 
profiles are given in Fig. 5, where the value of the phase 
varies continuously from — tt to tt, and the end point of 
the boundary between a tt phase line and a — tt phase line 
represents a phase defect (anticlockwise rotation denotes 
a vortex while clockwise rotation denotes an antivortex). 

For the cases of fully separated phases, we can see 
that there is a visible density hole in the density pro- 
file of the interior component [Figs. 4(a2) and 4(dl)] 
and a phase defect at the center of the corresponding 
phase profile [Figs. 5(a2) and 5(dl)]. The visible den- 
sity hole is referred to as a visible vortex (anticlockwise 
rotation) ^ because the topological defect is visible in 
the in situ density profile and it contributes to the an- 
gular momentum and the energy of the system. At the 
same time, there are three phase singularities that locate 



on the outskirts of the interior component and form a 
triangular lattice. Since these phase defects are invisible 
in the in situ density profile and contribute to neither 
the angular momentum nor the energy of the system, 
they are known as ghost vortices (exactly speaking, they 
are ghost antivortices because of their clockwise rotation) 
[37-39]. Of particular interest is the large density hole 
in the outer component [Figs. 4(al) and 4(d2)]. From 
its phase distribution [Figs. 5(al) and 5(d2)], we can 
see there are four phase defects in the region of density 
hole. The central phase defect is a vortex (anticlockwise 
rotation) and the peripheral three phase defects are an- 
tivortices (clockwise rotation) constituting a triangular 
lattice. The four topological defects are referred to as 
hidden vortices [37, 39, 40] because they carry significant 
angular momentum, though they are invisible in the in 
situ density distribution. Only after including the hid- 
den vortices can the well-known Feynman rule [s^, |43| 
be satisfied. The above analysis is also applicable to the 
cases of other states such as the inlaid separated-phase 
states and the mixed-phase states. For instance, for the 
partially mixed phases [Figs. 4(e), 4(f), 5(e), and 5(f)], 
there are four visible vortices in the density profile of 
each component in which one is a vortex lying on the 
trap center and the other are three antivortices locating 
on the periphery in a shape of triangular lattice. 

In general, the phase defects distribute regularly. The 
trap center is always occupied by a vortex in both species 
1 and species 2. For the partially mixed phases and the 
fully mixed phase, the outer antivortices in both compo- 
nents arrange in a triangular lattice at the almost same 
azimuth despite of different radiuses [Figs. 5(el)-5(f2)]. 
For the separated phases, the outer three antivortices 
in species 1 and those in species 2 arrange successively 
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along the azimuth direction due to the strong intra- and 
interspecies repulsions as shown in Figs. 5(al)-5(d2). 

As for the other cases of hk = hk 7^ hp = hp, our sim- 
ulation shows that there exist similar phase structures. 
To understand the formation mechanism of the above 
phase structures, we consider a single-component BEC 
with a VAVSS with arbitrary unequal winding numbers 
Ik 7^ Ip- As a general rule, we find that the steady state 
has N = lk-\-lp peaks and + 1 singularities in the phase 
distribution. When h > Ip, there are N singularities with 
unit positive topological charge on the periphery and Ip 
singularities with unit negative charge at the center, so 
that the total topological charge is s = N — Ip = h. 
When Ik < Ip, there are N singularities with unit neg- 
ative charge on the periphery and h singularities with 
unit positive charge at the center, so that the total topo- 
logical charge is s = N — h = Ip. For the case of h = ^p 
and = = 1/2, there is no phase singularity in the 
phase distribution due to the zero velocity field [28.]. If 
Ik = Ip and > (or a'^ < f3'^), there would be h 
single-quantum vortices (or antivortices) at the trap cen- 
ter as a result of the larger vortex ( ant i vortex) ratio a 

(/?)• 

For a nonzero-temperature BEC in the presence of 
thermal atoms, a vortex state is generally unstable due 
to the incoherent interactions between the BEC and the 
thermal atoms [4l|. For a zero-temperature BEC, a 
singly quantized vortex with repulsive interparticle inter- 
action is stable while the stability of a multiply quantized 
vortex is determined by the property and strength of the 
interaction [42, 43]. Here an interesting question is what 
state a VAVSS with h = 1, Ip = 2, and = = 1/2 
will develop into. Counter intuitively, we find that the ac- 
tual state is a vortex- ant ivortex cluster state consisting 
of a singly quantized vortex and three singly quantized 
antivortices. The density profile and the corresponding 
phase profile are similar to Fig. 4(e2) and Fig. 5(e2), re- 
spectively. The underlying physics is that in the presence 
of a single-quantum vortex the two-quantum antivortex 
only decays into three instead of two sing ly q uantized an- 
tivortices because the triangular lattice [43[ has the 
lowest energy and is the most stable. It is well known 
that the angular momentum of a single vortex depends 
on its winding number while the average angular mo- 
mentum per atom of a condensate is determined by not 
only the winding number but also the position of the 
vortex in the condensate. Combining with the conserva- 
tion requirement of the average angular momentum per 
atom, then the vortex and the three antivortices rear- 
range themselves to a stable spacial structure with the 
lowest energy in which the vortex locates on the trap cen- 
ter and the three antivortices distribute on the outskirts 
by means of a triangular lattice form. The similar anal- 
ysis can be generalized to the other cases. In Ref. (42^ , 
it is pointed out that a doubly quantized vortex may be 
stable in a BEC for certain regions of the interparticle 
interaction strength. In the present system, we do not 
observe similar doubly quantized antivortex in the equi- 



librium structures (Fig. 4 and Fig. 5) even we change the 
values of the intra- and interspecies interaction strengths 
in a large scope. The dynamic processes and details of 
how the VAVSS decay into the vortex- antivortex cluster 
states require further investigation. Incidentally, when 
hk = hp 7^ hk = hp our simulation shows that the 
equilibrium density profiles may suffer from distortion 
even deletion of petals especially for those in a separated 
phase due to the symmetry breaking of two group wind- 
ing numbers and the intra- and interspecies repulsions, 
which implies that the phase structures in this case will 
become quite complex. Further work would be necessary 
in order to understand the more complex structures in 
this case. 

The petal structure of a VAVSS termed initially in Ref. 
[24] remind us of an azimuthon which is proposed lately 
in nonlinear optics [45] and extended recently to BECs 
[4^. We show that there exist evident difference and 
certain relation between the VAVSS and the azimuthon. 
The difference is that our theoretical model and the rele- 
vant studies [23429} including the original literature [i^l 
on the VAVSS work in the laboratory frame while the 
investigations concerning the azimuthon work in the ro- 
tation frame |46| . Consequently, for the latter case 
there is an additional term with respect to angular veloc- 
ity in the stationary equation of system [46] . The relation 
is that the azimuthon is a special example of the VAVSS 
in view of their solution expressions. For the same wind- 
ing numbers of vortex and antivortex, the wave function 
of the VAVSS is given by (p{r, 0) - ae'^^ + pe'^e''^^ [see 
also Eq. (5) in [28]]. If we take the relative phase ^ = 
and introduce a parameter p = (a — /3)/(a + we can 
rewrite the wave function as ip{r^ 0) ~ cosilO) ^ ip sinilO) . 
Here the coefficient a + /3 has been absorbed in the nor- 
malization constant of the wave function, and the param- 
eter p should satisfy the two conditions of < |p| < 1 and 
+ = 1. In this context, the azimuthon expressed 
by ^ ~ cos{W) + ipsin{W) in Refs. [ii, |46[ is indeed 
a particular case of the VAVSS. Obviously, a VAVSS is 
not equivalent to an azimuthon when S 0. The sim- 
ilar analysis is also applicable to the cases of different 
winding numbers of vortex and antivortex. 



V. CONCLUSION 

We have numerically studied the exact 2D steady 
structures of two-component BECs with respective 
VAVSS. Depending on the winding numbers of vortex 
and antivortex and the intra- and intercomponent in- 
teractions, the system shows rich phase configurations 
such as fully separated phases, inlaid separated phases, 
asymmetric separated phase, and partially mixed phases 
with (deformed) petal-like component density distribu- 
tions. For given parameters, we display a phase dia- 
gram of two-components with respective VAVSS with 
hk = hp = ^ and hk = hp = 1- We show that the 
kinetic energy plays a key role in determining the struc- 
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ture of the two-component BECs with respective VAVSS, 
where the conventional phase separation criterion is in- 
apphcable and the TF approximation may fail due to 
the VAVSS, especially for the cases of VAVSS with large 
winding numbers. In addition, the conventional criterion 
for phase separation and the TF approximation can not 
discriminate different separated phases. According to the 
nonlinear stability analysis, the typical structures of the 
system are long-lived, which allows to be detected and 
tested in current experiments. The similar phase struc- 
tures also exist in the two-component BECs with respec- 
tive VAVSS with unequal winding numbers of vortex and 
antivortex. In this case, the density profile of each species 
forms a closed (highly modulated petal-like) structure ex- 
cept the fully mixed phase and the inner species in a par- 
tially mixed phase. In particular, an interesting vortex- 
antivortex cluster state occurs in each component in any 
of the possible phase structures, where the vortices and 



antivortices appear in the form of visible vortex, or hid- 
den vortex, or ghost vortex. Furthermore, a general rela- 
tion between the vortex-antivortex cluster states and the 
winding numbers of vortex and antivortex is revealed and 
analyzed. Finally, we show that the azimuthon proposed 
recently in nonlinear optics and in BECs is a special ex- 
ample of VAVSS in terms of their solution expressions. 
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